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Abstract 

The Mitsui model for the order-disorder type ferroelectrics is studied within 
the mean field approximation. The phase diagram of the model is obtained, 
its dependence on the transverse field is obtained. A scheme of setting the 
values of the model parameters is proposed; all the dielectric characteristics 
of Rochelle salt are calculated with the found parameters. Within the Mitsui 
model we study the physical characteristics of RbHS04 type crystals; the found 
values of the model parameters provide a good description of both dielectric 
and thermodynamic properties of these compounds. 
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1 Introduction 



The Mitsui model proposed [Q for description of Rochelle salt ||^ turned 
out to be suitable for description of the dielectric properties of other 
ferroelectric crystals (for instance of RbHS04 js) type) as well. This 
model is based on the assumption that the ordering structure elements 
move in double-minima potential wells. Due to the existing in these 
ferroelectrics internal electric field, directed to the opposite sides at the 
neighboring sites, this potential is asymmetric. In the case of Rochelle 
salt, spontaneous polarization is most likely caused by the hydroxyle 
group (0H)5 whereas in NH4HSO4 and RbHS04 spontaneous dipole 
moments are related to ordering of sulphate groups . 

For the first time, the pseudospin formalism for the Mitsui model 
was used in [Q. The values of the effective interaction parameters were 
found, which anebled a description of certain characteristics of Rochelle 
salt in the mean field approximation. It should be noted that the Mitsui 
model can have qualitatively different behavior depending on the values 
of the parameters; this dependence was explored in 1^. Infiuence of the 
model parameters on its thermodynamic properties was analysed; the 
phase diagram of the model was obtained; a region of possible values 
of the model parameters suitable for Rochelle salt was found. The di- 
electric properties of Rochelle salt were studied in 0, ||, || within the 
mean field approximation; however a description of the static dielectric 
susceptibility and the Curie- Weiss constant was not good. In |l^ the 
Mitsui model was treated within the two-particle cluster approximation; 
satisfactory results were obtained. On the other hand, in ||l2|, |l^ 
approximations of the interacting chains for the Mitsui model were pro- 
posed. 

The Mitsui model was also used for description of NH4HSO4 and 
RbHS04 crystals | 3|, |l5[ |l6|. Relaxation dynamics of these compounds 
was studied in 0, 17| |l8|, Unfortunately, within the mean field 

approximation it is impossible to obtain such values of the model param- 
eters which would predict the low-temperature first order paraelectric- 
ferroelectric phase transition in NH4HSO4 crystals It was shown 

that the temperature dependence of polarization in these crystals 
can be described with the temperature dependent model parameters. 
A more rigorous theory was developed in where exact taking into 
account of interaction between neighboring chains allowed to reproduce 
the correct temperature dependence of polarization in NH4HSO4. 

An important problem is to take into account tunneling of the or- 
dering elements through the potential barrier (a transverse field in the 
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Hamiltonian of the Mitsui model Q. This problem was studied in 
[|[ ^, ^ ^ . In Q the isotopic effect in deuterated Rochelle salt was 
described by changing tunneling only. In pl| the spontaneous polariza- 
tion and static dielectric susceptibility of the Mitsui model are calculated 
for several values of tunneling integral; a dynamic susceptibility from the 
Bloch type equation is obtained. 

It should be mentioned that in Rochelle salt a strong piezoelectric 
effect is present; hence for a correct description of all its characteristics 
we should take into account an interaction of the ordering structure 
elements with a phonon subsystem. The pseudospin-phonon interaction 



was taken into account in |24|, |25 



The aim of the present paper was to develop a fitting procedure for 
the ferroelectrics described by the Mitsui model and to find out whether 
their physical characteristics could be described within the mean field 
approximation. 



2 Thermodynamics of the Mitsui model in 
the mean field approximation 

Let us consider a Hamiltonian of ferroelectrics placed in external electric 
field and described by the Mitsui model |[ |j: 

^= E E J^,S^if)S■{f) - E ^..Sf (1)^1(2) 

-A5](Sf(l)-5f(2))-^i?5]E5r(/)-17EE5f(/), (2.1) 

i f i f i 

where Sf{f) = ^crf (/) is the a-component of the pseudospin operator, 
Jij and Kij are interactions between the spins of the same and of different 
sublattices; fx is the effective dipole moment of the ordering structure 
elements; A is the magnitude of asymmetry of the double-well potential; 
the transverse field characterises a finite height of the potential barrier. 

Neglecting quadratic in fluctuations terms, we obtain the correspond- 
ing Hamiltonian of the mean field approximation 



f 



(2.2) 



i?(l,2) = :^+Mi?T(^^-A). 



2 



Here £ = (^^(l)) + (5*^(2)) is the parameter of a ferroelectric ordering, 
whereas a = — {Sl{2y) is the parameter of an antiferroelectric 

ordering; J, K are the Fourier transforms of the interaction constants 
at q = 0. 

In order to calculate the partition function, we need to exclude the 
terms in the Hamiltonian containing the operators crf{f). With the help 
of the rotation transformation for the spin operators 



where sin^j = 
Hamiltonian 



^Kf) = '7^'(/)cOS</)/ +CTf(/)sin0/, 

^tif) = -'^r(/)sin<^/ +<(/) cos (t)f, 

Q 



COS d) f —— "^(■^) 



(2.3) 

we get the transformed 



H 



E 



J + K 



e+ 



J-K 



f 



Now we can easily obtain the free energy per site 



F 

N 



J + K 



J-K 



- lEln2coshM±^, 



(2.4) 



(2.5) 



where f3 = -^Kp is the inverse temperature. In the subsequent calcula- 
tions it is convenient to use the free energy normalized per "^^^ : 



/ 



F 



J + KN 



t 



In 2 cosh — + In 2 cosh — 



t 



.(2.6) 



Here we introduce the notations Ki 
V(C + a<^-7 + e)2+a;2, a = ^ 



v/(C-ocr + 7 + e)2+w2^ K2 



2A 
K+J' 



t = 



- K+J 

units. 



4T _ 2nE 



; the parameters with tildes, like J = J/ksi are in temperature 



Unknown parameters £, a are determined from the condition of ex- 
tremum of the free energy as a function of £, a. This yields the following 
system of equations 



CT + 7 + e K-i^ £ 

— tanh — H 

Ki t 



(7-7 + e K2 

—- tanh — 

K2 t 



— tanh — tanh — 

Ki t K2 t 



(2.7) 
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Here it is required that its solution (^o, '''o) is a saddle point of the func- 
tion /, that is, this function should have a maximum in the parameter 
of an antiferroelectric ordering o < 0) and a minimum in the pa- 
rameter of a ferroelectric ordering ^ < 0) [p6| . 

Using the known thermodynamic relations, we can calculate the other 
thermodynamic functions: entropy 

- i ^i^itanh^+A'atanh^^, (2.8) 
and internal energy 

u = , - — = / + ts = £2 _ aa^ -S^ KftSinh^. (2.9) 

j + KN ^ y t ^ ' 

To calculate the specific heat 

ds ( ds ds ds da 

" ^dt ^ ^ \dt ^ d^di ^ d^'dt 

we need to solve the system of equations for ^ and 
da fd^ da\ f d£ da\ 



where M, = ^ (iStzT + S ^^^^ ¥ j ' = i^l.'^L ■ Solution of 
this system reads 

da A^2 (Ml - 1) - Ni {M2 - 1) 



dt " {Ah - 1)(1 + aAh) + (M2 - 1)(1 + Ml) ' 

d^ _ N2{1 + aMi) + Nijl + aAh) 

dt ~ (A//i + aMa) + {A'h -!)(! + Mi) ' 

From this we obtain the specific heat: 

Kl K? ^ 
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cosh^ ^ cosh^ 



(1 - a)(iVi + + 2fl(M2iVf + MiTVl) ' 
(Ml ^ l)(aM2 + 1) + (M2 - l)(aMi + 1) 



(2.10) 
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Polarization of a crystal is calculated from the following thermodynamic 
relations: 

dF 

Experimentally, polarization per unit volume is measured: 

P 



P = 



V 



(2.11) 



where v is the unit cell volume. Hence, the static susceptibility is ob- 
tained in the following form: 



^ dE 



dp 



1 



V kB{K + J)de 

Af2(l+aMi) + Mi(l+aM2) 



(2.12) 



where ^ = 

"'^^^^^ de (l-Mi)(l+aM2) + (l-Af2)(l + aMi) • 

If we neglect tunneling {p. = 0), the expressions for the thermody- 
namic quantities and the system of equations for ^ and a simplify. The 
system of equations for the order parameter now reads 



<T = 



. K[ K' 

tanh h tanh — 

t t 



— - tanh 
t t 



tanh 



(2.13) 



where K[ = ^ — acr + 7 + e, K2 = ^ + acr — 7 + e. 

The expressions for the thermodynamic functions simplify as well. 
Entropy: 

s=ln ^2 cosh +ln ^2 cosh ~-^[2{^ + e)^-2{aa-^)a\, (2.14) 
internal energy: 



u = -{^ + 2e)^+{a(j-2-t)a. 



(2.15) 



3 Phase diagram of the Mitsui model 

It could be easily verified that the Hamiltonian and thermodynamics 
functions of the Mitsui model can be expressed via the dimensionless 
parameters a, 7, t, e, co. Depending on the values of these parameters, 
the model has different types of temperature behavior. 
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) the case J7 = {cu — 0). 

In order to explore the phase diagram in this case we use the scheme 
proposed by Vaks |^ . To find possible second order phase transitions in 
the system, we expand the system (2.13) in ^ up to the linear terms: 



^ = 1 - tanh 



2 -acr + j\ £, 



t J t' 

^^tanh^^^i^. (3.1) 

We obtain a system of equations for temperature of the second order 
phase transition. From the second equation we find that t = ^^fj^ and 
reduce (3.1) to a single equation for a: 

= l-a\ (3.2) 



-OCT + 7 -, 2 



arcthcr 

Let us introduce the notation a = tanh a and rewrite this equation as 

/(a) — —7 + a tanh a + a{l ~ tanh^ a) = 0, (3.3) 

a varies from to cx); the values of /(a) in the limiting points are 
/(a = 0) = —7, /(a — > 00) = —7 + a. Since 

d f 

/'(") = TT- = (a + l-2Q;tanha)(l-tanh^a) (3.4) 
da 

is a monotonic function, it turns to zero only at one point ao and has 
an extremum at this point. 

Let us consider the case 7 > a {f{a 00) < 0). The /(a) is an 
increasing function at small a and a decreasing one at large a. Hence, it 
has a maximum at certain ap. If f{ao) < 0, the model cannot undergo 
the second order phase transitions, whereas at /(ao) > there could 
be two phase transitions in the system. The line on the phase diagram 
(a, 7), separating these two characteristic types of behavior (with two 
phase transitions and with none) is determined from the following system 
of equations 

/(ao) = —7 + atanhao + q:o(1 — tanh^ ao) = 0, 

/'(ao) = (a + 1 - 2aotanhao)(l - tanh^ ao) = 0. (3.5) 

An equation for this line can be written in implicit form 



+ + ^-_±1, (3.6) 



6 



On increasing a, when a = 7 one of the second order phase transitions 
disappears, and the system goes to the region with a single second order 
phase transition. 

At T = the system can be in two states: the ordered = 1, ct = 0) 
or anti-ordered = 0, a = 1) one. Which of the two states is reahsed 
depends on the corresponding energies 

[/(^ = 1,CT = 0) = 1, U{^ = 0,a = l)=a-2-f. (3.7) 

In the region 

a - 27 > 1 (3.8) 

the system is in the ordered ground state since its energy is the lowest, 
whereas in the region 

a - 27 < 1 

the system is in the disordered state. 

The phase diagram of the model is depicted in fig.3.1. The line AB 
corresponds to the case a = 7, CED corresponds to equation (3.6), 
CEB corresponds to conditions (3.7) and (3.8). According to the ob- 
tained above results, in the region V of the model parameters values, 
the ferroelectric order in the system is absent at any temperature. In 
the region IV the ground state of the system in not ordered (^ = 0); 
however, here two phase transitions with the order parameter ^ are pos- 
sible. This means that ^ 7^ in a limited interval of temperatures, as 
shown in fig. 3. 2 for a = 0.5, 7 = 0.77. In the region II the system still 
can undergo two second order phase transitions, but its ground state is 
ordered. In this region several types of the system behavior are possible; 
they will be considered later. In the region I the system undergoes a 
single first order phase transition. We fail to obtain simple expressions 
for temperature of this transition and restrict our consideration by nu- 
merical analysis of temperature behavior of ^, in order to find new types 
of the temperature behavior of the model. The results of this analysis at 
different a, 7 (along the lines a = const at the phase diagram) are shown 
in figs. 3. 2-3. 8. We find that in a narrow layer Ila of the region II there 
emerges an additional first order phase transition at low temperature. 
The temperature dependence of ^ and other thermodynamic functions 
of the model in this region is shown in figs. 3. 4, 3.5 for the curves cor- 
responding to a = 0.3, 7 = 0.642. In the region lib the lower second 
order phase transition becomes of the first order. Qualitatively, the lines 
a = 0.3, 7 = 0.645 in figs. 3. 4, 3.5 illustrates the model behavior in this 
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Figure 3.1. (b) Phase diagram of the Mitsui model. 
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Figure 3.2: Temperature behavior of ^ at a = 0.5 and different values of 
7- 
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Figure 3.3: Temperature dependence of {d^/de) ^ at a = 0.5 and differ- 
ent values of 7. 
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Figure 3.4: Temperature dependence of ^ at a = 0.3 and different values 
of 7. 




Figure 3.5: Temperature dependence of (d^/de) ^ at a = 0.3 and differ- 
ent values of 7. 
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Figure 3.6: Temperature dependence of ^. 
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Figure 3.7: Temperature dependence of ^ at a = —0.3 and different 
values of 7. 
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Figure 3.8: Temperature dependence of (d^/de) ^ at a = —0.3 and 
different values of 7. 



regions. Deeper in the region II, the first order phase transition disap- 
pears, because the corresponding transition temperatures become equal 
(a = 0.3, 7 ~ 0.652 in figs. 3. 4, 3.5). On going from the region III to the 
region II, the first order phase transition persists up to the boimdary of 
the region 11c. The jump of the order parameter at this boundary turns 
to zero. At all other points of the phase diagram, the system undergoes 
a single second order phase transition. 

b) the case f2 ^ 0. As in the case O = 0, in order to find possible 
second order phase transitions at f2 7^ we expand the system (2.7) in ^ 
up to the linear terms. After some transformations, it can be presented 
as 



t = ^, i^^tanh.V(^^^-)^ 



[(7-a(7)2+a;2]3/2 t 
(7 - aay + y t j 
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Excluding t, we obtain the system (3.9) in the form f{a) = 0, where 



arcth[crvT" 



(7 - aa) - 



~ 2 
uj a 



1+U) 



r,2' 



(3.10) 



Here a) = a;/(7 — aa), a changes from up to l/\/l + w^. Let us find 
/(cr) in the limiting points 

/(c7 = 0) = -V7^T^< 0, 

1 -7 + ac7 + cr(l - cr^) 



(3.11) 



Hence, the line delimiting the region with a single possible second order 
phase transition is a solution of the equation f{amax) = 0. We obtain 
the following equation for o, 7: 

- 7 + ao- + (7(1 - CT^) = 0; (3.12) 

we also should solve the equation for a: a = + (D^. An equation 

for the line delimiting the regions with two possible second order phase 
transitions is obtained from the system of two equations = 0, 

/(cr) = 0. The expression for the derivative is too cumbersome to 

be presented here. The equations arc solved numerically; the calculation 
results are shown in fig. 3. 9. The line separating regions with ferroelectric 
ground state and a disordered ground state is found from the condition 



^ 0) = = 0). 



(3.13) 



This equation is also solved numerically. The phase diagram of the 

quantum Mitsui model at the value of tunneling oj = 0.4 is shown in 
fig. 3. 9 by dashed lines. From this figure one can see that the boundary 
of existence of the ordered ground state in raised up and now goes along 
the line CE'D'. The line of the two phase phase transitions corresponds 
to CE'B', whereas the line of a single phase transition corresponds to 
A'B'. 



4 Thermodynamics of Rochelle salt and 
deuterated Rochelle salt. 

The ordered phase in the Rochelle salt type crystals exists in a certain 
temperature interval: from Ti = 255K to T2 = 297K for a regular 
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Figure 3.9: The phase diagram of the quantum Mitsui model (solid and dashed lines correspond to 
uj = and uj = 0.4, respectively.) 



Rochelle salt and from Ti = 251K to T2 = 307K for deuterated Rochelle 
salt. The transitions to the ordered phase are of the second order in 
both crystals. According to the phase diagram (fig. 3.1) the possible 
parameters for Rochelle salt should be taken from the region IV [|[ ^• 
A difference between the transition temperatures Srs = {T2 — T'i)/(T2 + 
Ti) = 0.076 is small; therefore, the feasible values of the parameters 
for Rochelle salt are close to the line ED. In fig. 4.1 the line which 
corresponds to the sets of the theory parameters for Rochelle salt is 
depicted; above it there lies a line for deuterated Rochelle salt at which 
SdRs = 0.1. 

The parameters are chosen as follows: for the chosen at the phase 
diagram values of a, 7, we find such K, J, A for which the theoretical 
transition temperatures coincide with the experimental values; the effec- 
tive dipole moment /x is chosen such as the theoretical and experimental 
maximal values of polarization coincide. In figs. 4. 2-4. 5 we plot the exam- 
ples of calculated polarization and inverse static dielectric susceptibility 
for Rochelle salt and deuterated Rochelle salt at several different sets 
of the parameters. One can see, that the best fit for the susceptibil- 
ity is obtained at the small values of a and 7 from the curves Rs and 
dRs (fig. 4.1). Hence, the set of the model parameters which provides the 
best description of temperature behavior of Rochelle salt is: a = 0.29506, 
7 = 0.648 {K = 1473.59 K, J = 802.12 K, A = 737.33 K). The unit 
cell volume v — 5.24 x 10~^^ cm'^ of the model is twice smaller than 
the real one, sincethe full unit cell contains four ordering elements p7[ |, 
whereas the Mitsui model corresponds only to two ordering elements 
per unit cell. The effective dipole moment /x for Rochelle salt is found 
from the experimental value of polarization is ^ = 3.04 x 10^^^. For 
the deuterated Rochelle salt the best set of the fitting parameters is: 
a = 0.29952, 7 = 0.65 {K = 1502.83 K, J = 8102.07 K, A = 751.69 K, 
/z = 3.6 X 10^18). 

It should be noted that the Rochelle salt was studied in the mean 
field approximation by other authors too l^], ^ |[ 1). The parameters 
values for Rs and dRs obtained in these papers are given in Table 4.1. 

In order to explore the character of isotopic effect in Rochelle salt, 
we studied the influence of transverse field on the thermodynamic prop- 
erties of the model. Another fitting procedure for the Rochelle salt was 
to include tunneling into the model with the chosen already parame- 
ters for deuterated Rochelle salt, since deuterated mainly decreases the 
tunneling of the ordering structure elements. The calculation results are 
shown in figs. 4.6, 4.7. They indicate that it is impossible to describe the 
present isotopic effect by changing the transverse field only. A thorough 
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Figure 4.1: The phase diagram of the Mitsui model in region of possible values of the parameters for 
Rochelle salt. 



Table 4.1: The sets of the parameters for regular and deuterated Rochelle 
salt obtained in the mean field approximation (the sets 4 and 5 are found 
in the present work). 



No 


compound 


K, K 


J, K 


A, K 


a 


7 


1. 
2. 
3. 
4. 
5. 


Rs M 
dRs M 
dRs § 
Rs 
dRs 


1599.84 

1563.84 

1480 

1473.59 

1502.83 


764.64 

790.85 

800 

802.12 

810.07 


815.9 

789.7 

740 

737.3 

751.7 


0.353 
0.328 
0.298 
0.295 
0.299 


0.69 

0.671 

0.649 

0.648 

0.65 




Figure 4.2: Temperature dependence of polarization of Rochelle salt at 
different sets of the model parameters (■ are experimental points |^). 



18 




~l ' I ' I ' I ' I ' I 
220 240 260 280 300 320 



Figure 4.3: Temperature dependence of inverse static dielectric suscep- 
tibility of Rochelle salt at different sets of the model parameters (■ are 
experimental points |^). 




Figure 4.4: Temperature dependence of polarization of deuterated 
Rochelle salt at different sets of the model parameters (■ are experi- 
mental points 1^). 
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Figure 4.5: Temperature dependence of inverse static dielectric suscep- 
tibility of Rochelle salt at different sets of the model parameters (■ are 
experimental points pq|). 



numerical study of transverse field effects on the thermodynamic char- 
acteristics of the model has shown that taking into account of tunneling 
does not improve an agreement between the theory and experiment. For 
a fixed 7, an increase in w only increases the static susceptibility. 

The performed calculations of spontaneous polarization and static 
dielectric susceptibility for Rs and dRs have shown that it is impossible 
to find such sets of the microparameters K , J, A, and /i, which would 
provide a simultaneous good description of both characteristics. There- 
fore, the statements of other authors that the Mitsui model is adequate 
to Rs or dRs and properly describes the physical characteristics of these 
crystals are incorrect. Within the model (2.1) either with tunneling or 
without it, this is not possible. 



5 Relaxation dynamics of Rs and dRs. 

It has been mentioned in several papers ^ that as indicated by 
measurements of dielectric dispersion, tunneling effects do not play an 
important role in Rochelle salt crystals. Therefore, the time correlations 
are believed to be caused by interaction with a phonon subsystem. A 
phenomenologic description of this interaction is performed within the 
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Figure 4.6: Temperature dependence of polarization of Rochelle salt and 
deuterated Rochelle salt for different sets of the parameters: lines for Rs 
and dRs correspond to sets of parameters 4, 5 in table 4.1, dashed line 
is the model with the set of parameters 5 and fl = 50. 
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Figure 4.7: Temperature dependence of inverse static susceptibility of 
Rochelle salt and deuterated Rochelle salt for the same sets of the pa- 
rameters as in previous figure. 
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Glauber method J32]. A master equation for a in this approach can 
be presented as g: 



d 1 
—a— (J = cr — - 
2 



tanh h tanh , 



t 



tanh — - — tanh — =■ , 
t t 



(5.1) 



where a is the parameter of the spin-phonon relaxation. 

Considering a weakly non-equilibrium system, we expand ^, a near 
their equilibrium values ^, ct: ^ = ^ + ^(i), ct = + cr(t). From (5.1) we 
obtain that i^, a satisfy equation (2.13), whereas £,{t), a{t) can be found 
from the following system of equations 



a^jT^it) = Mii^it) + Mi2a{t) - Mie(t), 
dt 

-a^a{t) = M2ii{t) + M22(T{t) - M2e(t), 
dt 



(5.2) 



where Mn = 1 - i(l - - f^'), M12 = "X^^' ^^21 = M22 = 

1 + £(1 „ ^2 „ ^2)^ ^ „ ^2 „ ^2)^ ^ _2i._ 

In a frequency representation, a linear response of ^(i) to a small field 
e{t) is obtained in the following form 

?(^) = Xo(^)e(w), 



Xo(t^)= 



WlTl W2T2 



1 + Lj2t-2 1 + t^2^2 
Mi{Mil - + M2M12 



^^2T-2" 

- Cj2t2 ^ 1 + Cj2^2 



(5.3) 
(5.4) 



W2 



Mi{Mil - + M2M12 
a^irr'-r2') ' 



_1 _ A/i 1 +M22 ib.^ (Ml 1 - Af 22 )^ +4Mi2 M21 



are the inverse relaxation 



where 2 - 2a 
times of the model. Relations (5.3), (5.4) for polarization P and electric 
field E now read 



2^2 



1 



:Xo(w). 



(5.5) 



V {K + jy 

The dielectric permittivity is related to the corresponding susceptibility 



as 



e{uj) = 1 + inxicj) - s\cj) - ie"(tj). 



(5.6) 
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Figure 5.1: Temperature dependence of e' for Rochelle salt; lines are 
theoretical results for a = 0.9 x 10~^^, ^ ~ 1.8487 x 10~^^; symbols are 
experimental points |^0|| . 

where e'(w) £"(w) are the real and imaginary parts of the dielectric 
permittivity. 

In figs. 5. 1-5.6 we depicted the calculated dielectric permittivity of 
Rochelle salt. In calculations we use the parameters found above: 
K = 1473, J = 802.12, A = 737.33 for Rochelle salt and K = 1502, 
J = 810.07, A — 751.69 for deuterated Rochelle salt. Unfortunately, at 
the already set values of fi we are unable to find a value of the relaxation 
parameter a such as to obtain a satisfactory description of the experi- 
mental points for the dynamic permittivity. Therefore, in this case we 
take /i = 1.8427 x IQ-^^ for Rochelle salt and ^ = 2.0 x 10"^** for deuter- 
ated Rochelle salt. Polarization and static susceptibility calculated with 
the new values of fi are shown in figs. 5. 7, 5.8. 

As one can see in figs. 5.1-5.6, the mean field approximation yields 
a fair description of the temperature behavior of the dynamic dielectric 
permittivity of Rochelle salt at all temperatures except for the vicinity 
of the phase transitions, where the real part of the permittivity tends to 
EoQ. This is due to the divergence of the relaxation times of the model 
at the transition temperatures (fig. 5. 9). 
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Figure 5.2: Temperature dependence of e" for Rochelle salt. 

6 Thermodynamics and relaxational dy- 
namics of RbHS04. 

Ferroelectric properties of RbHS04 were discovered in 1960 |4^. This 
crystal undergoes a single second order phase transition (Tc = 265 K). 
A microscopic theory of RbHS04 for the first time was proposed in 
|. The studied model was based on the assumption that the phase 
transition in RbHS04 is related to ordering of sulphate groups, which 
in a high-temperature phase move in asymmetric double- well potentials. 
Hence, RbHS04 can be described within the Mitsui model (see (2.1)). 
In 1^, |l^, |l^, |l^, |l^, the thermodynamic and dielectric properties of 
RbHS04 RbDS04 and disordered mixtures RbHi„a;Da;S04 were studied. 
Unfortunately, the found in these papers values of the microparameters 
do not provide a proper description of the experimentally observed jump 
of the specific heat at the transition point (Ac = 9.02 J/Mol K). 

Since RbHS04 and related disordered mixtures undergo a single sec- 
ond order phase transition, the possible values of the model parameters 
for these compounds should be taken from the regions I, II of the phase 
diagram. Also, since the specific heat of the studied model depends only 
on the dimensionless parameters a and 7, the values of these parameters 
were chosen such as to fit the theoretical jump of specific heat to the ex- 
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Figure 5.3: The frequency dependence of the real part of dynamic di- 
electric permittivity e'(i^) at different temperatures T (K): a) - 235, b) 
- 245,_c) - 265, d) - 285, e) - 305, f) - 315. Solid lines are the results 



from [M 
are taken from 



dashed line is the result of present paper, experimental points 
|3|gl,T-||, + -@,.-[g, x-||, A 



- |30 
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Figure 5.4: The frequency dependence of the imaginary part of dynamic 
dielectric permittivity e"{v) at different temperatures T (K): a) ~ 235, 
b) - 245^ c) - 265, d) - 285, e) - 305, f) - 315. SoHd Hues are the results 
from p9[ , dashed line is the result of present paper, experimental points 
are taken from ■ - O- lEsI, Ri, ▼ - [H, + - il, o - BtI, x - M, □ 
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Figure 5.5: Temperature dependence of e' for deuterated Rochelle salt; 
lines are theoretical results for a = 0.8 x 10"^^, /x = 2.00457 x 10"^^; 
symbols are experimental points 




Figure 5.6: Temperature dependence of e" for deuterated Rochelle salt. 
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Figure 5.7: Temperature dependence of polarization of regular (solid 
line) and deuterated (dashed line) Rochelle salt for sets of parameters 4, 
5 of table 4.1; = 1-8487 x 10-^^ HdRs = 2.0046 x 10~^^. 
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Figure 5.8: Temperature dependence of inverse static dielectric suscep- 
tibility of regular (solid line) and deuterated (dashed line) Rochelle salt 
for the same sets of parameters as in previous figure. 
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Figure 5.9: Temperature dependence of inverse relaxation times of reg- 
ular (solid line) and deuterated (dashed line) Rochelle salt. Symbols are 
experimental points of l3l|]). 



perimental value. Then the parameters if, J, A were determined from 
a and 7 using the condition that the calculated transition temperature 
should coincide with the experimental one. Finally, for RbHS04 we ob- 
tained the following set of the parameters K ~ 880.64 K, J = 780.9 K, 
A = 398.8 K, n = 0.425276 x lO'^^. Here we took into account the fact 
that for RbHS04 v = 2.105 x IQ-^^ ,,^^3 

Thermodynamic and dielectric characteristics of the disordered com- 
pounds RbHi_j;D2:S04 were calculated in the mean field approximation 
under assumption that the corresponding model parameters are given 

by 

Jx = Jh{1 - x) + Jdx, 
A^ = Aff(l-x)-HAi5x. (6.1) 

Here Kh, Jh, and K]j, Jjj, Ao are the model parameters, cor- 
responding to RbHS04 and RbDS04, respectively. We took into ac- 
count the fact that for RbDS04 Tc = 252 K. Analogously to the 
case of RbHS04, we obtained the values of the model parameters for 
RbH0.3D0.7SO4, and then using (6.1) found the parameters for RbDS04. 
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Table 6.1: The sets of the model parameters for RbHS04 and RbDS04, 
obtained within the mean field approximation (the sets 5 and 6 are found 
in the present work). 



No 


compound 


K, K 


J, K 


A, K 


a 


7 


1. 


RbHS04 




616 


784 


245 


-0.12 


0.35 


2. 


RbHS04 i 


16 


1 


1040 


228 


468 


0.64 


0.738 


3. 


RbHS04 i 


17 


1 


616 


777.2 


244 


-0.116 


0.35 


4. 


RbHS04 


15 


1 


616 


784 


244 


-0.12 


0.348 


5. 


RbHS04 






880.64 


780.9 


398.8 


0.6 


0.48 


6. 


RbDS04 






842.4 


747.1 


381.5 


0.6 


0.48 



Table 6.1 contains values of the model parameters for RbHS04 and 
RbDS04, obtained in the present and previous [|[ |l6|; |l^, ^ works. 
Using equation (6.1) and the presented in Table 6.1 values of the model 
parameters for RbHS04 and RbDS04, we calculated the physical charac- 
teristics of the mixed compounds RbHi_2:D2;S04. The obtained results 
for polarization, inverse dielectric susceptibility, and contribution of the 
ordering structure elements to the specific heat of RbHi_2;D2,S04 along 
with the available experimental data are shown in figs. 6.1-6.3. The 
fig. 6. 4 illustrates the calculated specific heat and the experimental points 
for RbHS04. In figs.6.5-6.9 we plotted the calculated s'{iy,T), e"{v,T) 
for RbHS04 and the inverse relaxation time RbHi_j.D2,S04. A good 
quantitative fit of the obtained theoretical results to experimental data 
for the studied characteristics of RbHS04 is obtained. If figs. 6. 10-6. 13 we 
show that the proposed theory also well describes the experimental data 
of Q for temperature and frequency dependences of e'{iy,T), e"{v,T) 
for RbH0.3D0.7SO4. Hence, a good agreement between the theoretical 
results and experimental data is obtained for thermodynamic and dy- 
namic characteristics of the RbHS04 RbH0.3D0.7SO4 ferroelectrics. We 
can maintain that the Mitsui model (2.1) is totally adequate to these 
crystals. 

7 Concluding remarks 

In this paper we calculated thermodynamic and dynamic characteristics 
of ferroelectric order-disorder type compounds described by the pseu- 
dospin model with asymmetric double-well potential. For the first time, 
a role of tunneling on the systems described by this model is studied 
thoroughly; its influence on the phase diagram is shown. Possible phase 
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160 180 200 220 240 260 280 

Figure 6.1: Temperature dependence of polarization of RbH3^Di_2:S04 
for X =1, 0.7, 0.5, 0.25, 0; □ are the experimental points taken from pTj . 
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Figure 6.2: Temperature dependence of the inverse static susceptibility of 
RbHa;Di_a;S04 for x =1, 0.7, 0.5, 0.25, 0; Symbols are the experimental 
points taken from |^ . 
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Figure 6.3: Contribution of the ordering elements to the specific heat as 
a function of temperature for RbH2;Di_2:S04 at x —\, 0.7, 0.5, 0.25, 0. 
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Figure 6.4: Specific heat of RbHS04 as a function of temperature (the 
sohd hue represents theoretical results obtained within the mean field 
approximation, the symbols are experimental points 0). 
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Figure 6.6: Temperature dependence of e" of RbHS04 {a = 0.48 x 
10"^'^): the solid line represents theoretical results; the symbols are ex- 
perimental points p^ . 
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Figure 6.7: e' of RbHS04 as a function of frequency; the solid lines 
represent theoretical results; the symbols are experimental points [ p7| . 
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Figure 6.8: e" of RbHS04 as a function of frequency; the solid Unes 
represent theoretical results; the symbols are experimental points [ p7| . 
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Figure 6.10: Temperature dependence of e' for RbH0.3D0.7SO4 [a = 
0.54 X 10~^^): the solid lines represent theoretical results; the symbols 
are experimental points . 
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Figure 6.11: Temperature dependence of e" for RbH0.3D0.7SO4 (a — 
0.54 X 10"^^): the solid lines represent theoretical results; the symbols 
are experimental points | pT[ . 
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Figure 6.12: e' of RbH0.3D0.7SO4 as a function of frequency; the solid 
lines represent theoretical results; the symbols are experimental points 
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Figure 6.13: e" RbH0.3D0.7SO4 as a function of frequency; the solid lines 
represent theoretical results; the symbols are experimental points 

transitions are analysed at different values of the model parameters; 
at these values the physical characteristics of the model are calculated. 
Interesting results are obtained for compounds undergoing three phase 
transitions, like {(NH4)3H(S04)2}i-:e{(ND4)3D(S04)2}:. !§■ Because 
of a limited number of experimental data available for these crystals, we 
cannot draw final conclusions about the model adequacy for these crys- 
tals. We also calculated the physical characteristics of the {dRs)xRsi-x 
and RbHi_2.D2,S04 type ferroelectrics. 

Our study has shown that within the studied model it is most likely 
impossible to describe the phase transition and physical characteristics 
of NH4HSO4 crystal. Interesting results are obtained for the RbIIS04 
and Rs type ferroelectrics. The first crystal is not piezoelectric in the 
paraelectric phase, whereas the piezoelectric effect determines the phys- 
ical characteristics of the second one The considered here model 
does not take into account the piezoelectric interaction. 

It should be also noted that in the mentioned ferroelectrics, an in- 
teraction of the ordering structure elements with phonons is rather im- 
portant. Corresponding studies for these crystals were performed in 



Hence, despite a great attention to the crystals described by the 
Mitsui model, especially to Rochelle salt and of RbHS04 type, it still 
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unclear whether this model provides a quantitative description of the 
entire spectrum of physical characteristics for these crystals. Indeed, in 
practically all earlier paper, the microparameters were chosen by fitting 
to only a few characteristics of the crystals, whereas other quantities 
were not calculated with these parameters. Naturally, the incorrectly 
chosen values of the K, J, A parameters led to the incorrect values of /z 
and a as well. The aim of the present paper was to calculate all physical 
characteristics of the Mitsui model, study their dependences on the mi- 
croparameters, explore the tunneling effects on these characteristics. We 
also intended to develop a fitting procedure for regular and deuterated 
Rochelle salt and RbHS04 type ferroelectrics, which would provide a 
good fit to the available experimental data for these crystals and help to 
establish the adequacy (or inadequacy) of the considered model to the 
specific crystals. 

In it has been shown that an important characteristics in the fit- 
ting procedure for ferroelectric materials is a specific heat. In the present 
paper at setting the microparameters K, J, A we used the data for the 
specific heat of the studied crystals. However, it should be mentioned 
that for specific heat Rs and dRs reliable data are most likely absent. 

The performed calculations yielded the values of the microparam- 
eters for the RbHS04 crystals, with which we calculated the physical 
characteristics of the disordered ferroelectrics RbHi_2;Dj;S04. The pro- 
posed consistent fitting scheme provides a good quantitative description 
of experimental data for pure crystals; we also present the results for dis- 
ordered crystals RbHi_2;Da;S04 which can be experimentally verified. 

For Rs and dRs crystals it is impossible to find such values of the 
model parameters, which would provide a simultaneous fit to the experi- 
mental data for Ps{T) and e(0, T). An agreement could be somewhat im- 
proved by taking into account the piezoelectric interaction. This is also 
indicated by calculations of the dynamic characteristics of the crystals 
for which a satisfactory agreement with experimental data is obtained. 
It is clear that with increasing frequency, the effects of piezoelectric in- 
teraction in Rs type crystals become unimportant. The crystals become 
clamped, and we obtain a good agreement with experimental data for 
e'{v,T) and e"{v,T). 

Another important point is, if we choose the values of the micropa- 
rameters K , J, A at which (see the phase diagram) the third low- 
temperature phase transition in Rs emerges, an agreement with the 
experimental data simultaneously for Ps{T) and e{0,T) could be much 
improved. However, the question about the low-temperature transition 
in Rs is still open, even though such a transition was reported in some 
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papers. 
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